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ABSTRACT 


Context. Stellar activity induced by active structures such as stellar spots and faculae is known to strongly impact the radial velocity 
(RV) time series. It is thereby a strong limitation to the detection of small planetary RV signals, such as that of an Earth-mass planet 
in the habitable zone of a solar-like star. In a serie of previous papers, we studied the detectability of such planets around the Sun 
observed as a star in an edge-on configuration. For that purpose, we computed the RV, photometric and astrometric variations induced 
by solar magnetic activity, using all active structures observed over one entire cycle. 

Aims. Our goal is to perform similar studies on stars with different physical and geometrical properties. As a first step, we focus on 
Sun-like stars seen with various inclinations, and on estimating detection capabilities with forthcoming instruments. 

Methods. To do so, we first parameterize the solar active structures with the most realistic pattern so as to obtain results consistent 
with the observed ones. We simulate the growth, evolution and decay of solar magnetic features (spots, faculae and network), using 
parameters and empiric laws derived from solar observations and literature. We generate the corresponding structures lists over a full 
solar cycle. We then build the resulting spectra and deduce the RV and photometric variations, first in the case of a “Sun” seen edge-on 
and then with various inclinations. The produced RV signal takes into account the photometric contribution of spots and faculae as 
well as the attenuation of the convective blueshift in faculae. We then use these patterns to study solar-like stars with various inclina¬ 
tions. 

Results. The comparison between our simulated activity pattern and the observed one validates our model. We show that the inclina¬ 
tion of the stellar rotation axis has a significant impact on the photometric and RV time series. RV long-term amplitudes as well as 
short-term jitters are significantly reduced when going from edge-on to pole-on configurations. Assuming spin-orbit alignment, the 
optimal configuration for planet detection is an inclined star (i 45°). 

Key words. Stars : activity - Stars : starspots - Stars : solar-type - Techniques : radial velocities - Techniques : photometric 


1. Introduction 

Thanks to the radial velocity (RV) technique, more than 500 
exoplanets have been discovered in the two last decades 
(http://exoplanets.eu). First limited to Jupiter-like planets on 
short period orbits, the RV technique allows today to detect 
Neptune-like planets (10-40 MEanh) and Super Earths (1.2-10 
MEarth) on longer period orbits. This progress was made pos¬ 
sible thanks to the important improvements on the sensitivity 
and stability of instruments, as well as obser vational strategies 
to aver age out known sources of stellar signals (iDumusaue et al.L 
1201 lah . For example, the High-Accuracy Radial velocity Planet 
Searcher (HARPS) s pectrograph locate d on the ESO 3.6 m tele¬ 
scope in La Silla (IPene et ^ l2002h gives a precision of 1 
m s ' for a Solar-type star in average conditions, compared to 5- 
10 ms * with previous instruments. Furthermore, future instru¬ 
ments, such as ESPRESSO on the Very Large Telescope (VLT) or 
G-CLEF on the Giant Magellan Telesco pe (GMT) are expected 
to reach a p r ecisio n down to 0.1 ms * dMegevand et al.U2010l 
iFzresz et al.L l2014l) . Such a level of precision will theorically 
give access to lower planetary masses far from their host stars 
such as Earth-like planets in the habitable zone (hereafter HZ). 

However, low-amplitude RV planet signals (such as RV sig¬ 
natures of Earth-mass planets in the HZ) are much more sen¬ 
sitive to stellar perturbations than giant planet RV signatures. 


In the case of a low-mass planet RV signature with an ampli¬ 
tude in the 0.1-1 ms ' range, the stellar noise or “jitter” is high 
enough to either mask or mimic the planet-induced RV varia- 
tions, even in the case of a chromospherica l ly nearly inactive star 

(llsaacson & Fisch^l2010ilLagrange et al.Ll2010l.lMeunier et al.L 

I2010h . 

The so-called stellar jitter mainly comes from three different 
sources: stellar oscillations or pulsations, granulation, and stellar 
magnetic activity. 

Stellar pulsations or oscillations dominate the RV jitter on 
the shorter timescales. For FGK dwarfs, they are mainly driven 
by accoustic or pressure waves (p-modes). Pressure waves are 
commonly believed to originate from turbulent convective mo¬ 
tions occuring in the stellar outer layers and propagate through 
the star. They induce RV shifts with an amplitude of a few 
cm s ' to one ms“'. The oscillation periods range from a few 
minutes (e.g., five minutes in the case of the Sun) to a few tens 
of minutes. The amplitude and the period of the oscillations in¬ 
crease with stellar mass. Long exposure times and observational 
strategies have shown to be efficient in averaging the RV noise 
induced by pulsations in the ca se of late-type stars (ISantos et al.L 
|2004 IDumusaue et al.Ll201 Ibl) . 

The photospheric granulation phenomenon accounts for the 
convective plasma motions occuring in the outer envelop of 
solar-mass stars. The quick rise and fall of the plasma results 
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in a pattern of bright granules and darker lanes at the stellar 
surface, and in rapidly evolving RV shifts, with amplitudes up 
to a few km s ' locally. When integrated over the entire stel¬ 
lar disc, upflows and downflows average out, leaving a resid¬ 
ual RV jitter at the level of the ms *. Up to now, there have 
been fe w tentatives to e s timate t he impact of granulation on 
the RV. iDumusaue et al.l (l201 Ibl) made a first estimation for 
different spectral type stars based o n the study of asteroseis- 
mology measurements. ICegla et aP (12013h developped a four- 
component model of granulation, building absorption line pro¬ 
files from three-dimensional magneto hydrodynam i c sola r sim¬ 
ulations. This study was continued by ICegla et al.l (l2014l) . who 
derived the corresponding RV time series as well as other time 
series of observables such as the bisector amplitude. They con¬ 
cluded first that granulation has a strong impact on RV at the 
cm s * to m s * levels and would be a potentially significant lim¬ 
itation to low-mass planets detection. Second, they found corre¬ 
lations between the RV time series and other observables (e.g. bi¬ 
sector curvature or bisector inverse slope) that already allow 
to partially correct the granulation signal and reduce the corre¬ 
sponding jitter by one third. Meunier et al {submitted) simulated 
a collection of both granules and supergranules on the solar visi¬ 
ble hemisphere over a whole solar cycle and derived the induced 
RV and photometric time series. These authors concluded that 
the granulation noise is hard to average out even over an entire 
night and will therefore have a significant impact on detection 
limits. 

From a few days to several years timescales, the stellar jitter 
is dominated by the so-called stellar magnetic activity. Active 
regions such as dark spots an bright faculae (colder and hot¬ 
ter than the quiet photosphere, respectively) grow and evolve 
on the stellar surface, inducing changes in the stellar irradiance 
and RV variations (as the flux loss or excess in the active re¬ 
gions distorts the CCF of the stellar spectra). The stellar rota¬ 
tion makes one see these active regions moving across the stellar 
disc, thus inducing an apparent Doppler shift. The first studies 
of the impact of sta r spots on RV and l i ne bisectors were made 
by S^r & Donahu'3 (Il997l) and lHatz^ (l2002h . lSaar & Donahue! 
(119971) built a semi-empirical law that directly bound the spot- 
induced RV semi-amplitude, the spot fi lling factor and th e stel¬ 
lar rotational velocity for cool F-G stars. ISaM(l2003ll2009h made 
similar stu dies for bright facu lae. The first spot model was de¬ 
veloped by |Desortet^(|2003)- The authors computed synthetic 
stellar spectra and applied a black-body law to take into account 
the contribution of a colder spot. They quantified accurately the 
RV amplitude, RV bisector and photometric variations induced 
by a starspot for a range of different stellar and spot properties. 

Such simple spot models have proved to be useful and effi¬ 
cient to simulate the RV signature of one or a few active regions 
during a timescale of the order of the stellar rotational period. 
They are well adapted to stars that host a single main spot or a 
few ones, such as active young dwarf stars {e.g. HD 189733, see 
lDumusau3.120l4 and ref. therein). They are however not suffi¬ 
cient to reproduce the total activity-induced RV signal for Sun¬ 
like stars, due to the presence of multiple evolving active regions 
on the stellar disc, and d ue to the strong contribution of faculae 
becau se of slow rotation (iMeunier et al.ll201^ IDumusaue et al.L 
l2014l) . A complete model of the stellar activity pattern for longer 
timescales (of the order of the activity cycles) is needed. For 
most stars though, we have no detailed information on the ac¬ 
tivity properties as active structures are not directly observable. 
Doppler and Zeeman-Doppler imaging technics allow to recover 
the largest active structures (or active structure clusters) only 
for fast-rotating stars {i.e., young, active ones) but not for old. 


slow-rotating solar-type stars (lRiceLl2()()^ . The structure proper¬ 
ties, as well as the convection ones, are mostly unknown and we 
rely only on indirect and global estimators s uch as the Calcium 
(Ca) index (see e.g. Baliunas et al.L 119831 iNoves et akl 11984 
EaU & Lockwood! 1200^ 


The Sun is an exceptional star, for which we have a lot more 
of information. Indeed, both solar dark spots and bright features 
are well observed and the various solar activity properties are 
thus well described. The Sun, if seen as a moderately active star, 
therefore represents an ideal prototype to study the impact of ac¬ 
tivity on low-mass planets detectability. In a series of papers, we 
modelled with a great accuracy the solar activity pattern using 
detailed observations {i.e., spots catalogs and magnetograms) ex¬ 
tending over the full solar Cycle 23. We then rebuilt the induced 
photometric variations and the corresponding spectra, and de¬ 
duced from t he latter the RV t ime se ries. We first considered only 
cold spots in lLagrange et al.l(l2010t) (hereafter Paper I), studying 
their impact on the detec tability of an Earth-m ass planet in the 
HZ of a Solar-like star. In IMeunier et akl (1201 Ol) (hereafter Paper 
II), we extended this study by considering in addition the con¬ 
tribution of bright faculae and, for the RV, the attenuation of the 
convective blueshift in magnetically active regions. Combining 
the effect of the three activity components, we found out that the 
attenuation of convective blueshift in faculae (hereafter the con¬ 
vective component) dominated the induced RV variations (this 
happens as the Sun is a slow rotator: for larger v sin i, the pho- 
tospheric effect of spots becomes dominant). We showed that 
unless with correction tools (still to be identified and tested) it 
would be impossible to detect an Earth- like planet in the HZ . 
Using the same simulation as in Paper II. lLagrange et al.l (1201 ll) 
(hereafter Paper III) estima ted the astrometric effect of stellar 
magnetic activity. Finally, in iMeunier & Lagrange! (l2013l) (Paper 
IV), we used the strong correlation between the convective com¬ 
ponent and the Ca index (as both are directly related to the chro¬ 
mospheric plage or photospheric facula filling factor) as a tool 
to correct RV time series from the convection-induced jitter. We 
showed that an Earth-like planet in the HZ would become de¬ 
tectable given that: i) there is an excellent signal-to-noise ratio 
(hereafter S/N) on the Ca index data; and ii) the temporal cover¬ 
age of the RV observations over one activity cycle is very dense 
(typically one night out of four during the cycle). 


The objective of the present paper is to extend the previ¬ 
ous study to a Sun-like star seen under any inclination. To do 
so, we fully parameterize the activity pattern of a Sun-like star. 
A few activi t y para metrizations haye already been deyeloped. 
iBarnes et al.1 (l201 ill undertook the first tentatiye to parameter¬ 
ize dark spot distributions at low and high actiyity leyels in 
the case of M-type stars. A similar simple model was used by 
Lleffers et~^ (l2014l) . this time for young and actiye G and K 
dwarfs. A parametrizati on of an actiyity pattern was made by 
IDumusaue et akl (1201 lal) . who simulated dark spot distributions 
oyer the stellar surface for different leyels of actiyity of a solar¬ 
like star. The authors deriyed detection limits for different actiy¬ 
ity leyels and obseryational strategies. This study was howeyer 
limited to dark spots and did not include the flux effect of bright 
features (faculae and network), nor the inhibition of the conyec- 
tiye blueshift in actiye structures. 

Using the Sun as a template presents seyeral important ad- 
yantages: i) we haye much more information on solar actiyity 
than for any other star; if) the yalidity of such a model can be 
easily asserted by a comparison with the preyious papers; and 
Hi) such an actiyity model can be adapted to other spectral types 
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and other activity levels than the Sun, if the active region con- 
hguration is supposed similar. In addition, it allows studying the 
impact of stellar inclination on the results, contrary to the solar 
observations for which only one hemisphere is observed at any 
given time. 

In Sect. |2] we build the spot and facula distributions on the 
Sun surface over a full solar cycle. We describe in details the 
laws and parameters we used to make these simulations. To test 
our approach, we then compare the obtained structure distribu¬ 
tions with the observed solar distributions used in paper II. In 
Sect. 13 we compute the induced RV variations, taking into ac¬ 
count both the photometric contribution of spots and faculae and 
the attenuation of the convective blueshift in active structures. 
Then we compare our results with those previously obtained us¬ 
ing observed solar magnetic features. In Sect. |4] we study the 
influence of the stellar inclination on the resulting RV and pho¬ 
tometry. The astrometric time series will be treated in a separate 
paper. We also compute the corresponding detection limits. We 
hnally discuss our results and the use of our model for other stars 
in Sect.|5] 

2. Building the activity pattern 

2.1. Approach 

We consider the solar activity pattern during Cycle 23. We use 
the Sun as a template and as a mean of comparison with the 
observed activity pattern used in Papers I and II. We parameter¬ 
ize an extended range of activity levels and timescales. We start 
from the general distribution of the activity level over the cy¬ 
cle, then we parametrize the spatial and temporal distributions 
of active structures and their dynamics, and we hnally mod- 
elize the individual behaviour of the structures, including the 
dark spots, the bright faculae around these spots and the net¬ 
work smaller features. The two outputs of our simulations are 
the lists of the dark spots and of the bright features, respectively. 
Each hie gives, for each time step of the simulation, the structure 
sizes (in millionth of hemisphere, hereafter /iHem), latitudes and 
longitudes. For a given set of parameters, these hies will then 
represent the inputs of our simulation tool to produce the spectra 
and the corresponding RV at each time step, as it was done in 
Papers I and II with observed solar structure lists. 

2.2. Building the spot and facula catalogs: input parameters 

All input parameters are summarized in Table. [T] We detail them 
hereafter. 

1. Global activity level. The sunspot-induced solar activity fol¬ 
lows a cycle of 12 + 1 years. To mimic as accurately as pos¬ 
sible the global activity level and its distribution over such a 
cycle, we used the relative sunspot number, or Wolf number 
(noted R), as a proxy. For this study, we based our simula¬ 
tions on solar Cycle 23. This allowed us both to build a re¬ 
alistic activity pattern and to compare our results with those 
previously obtained in Paper II. The daily Wolf number R 
for Cycle 23 was recovered at the Solar Inhuences Analysis 
Data Centefl (SIDC) and is displayed in Fig.[T](black dots, 
upper panel). We hrst smoothed the daily R data to obtain a 
long-term reference Rsmth (Fig- El solid line, two top panels) 
for the cycle shape. 

The observed daily dispersion of the Wolf number comes 
from two contributions: the hrst comes from the rotational 

' available at: http://sidc.oma.be/ . 





Period (doy) 

Fig. 1. Top: daily {dots) and smoothed (bold line) Wolf number 
over solar cycle 23. Middle: same smoothed Wolf number (black 
line), and when randomly adding a second order polynomial dis¬ 
persion (red dots). Bottom: Power spectra of the observed daily 
(black) and simulated (red) Wolf numbers (respectively R and 
R ). The solar rotation period (blue solid line) and its hrst har¬ 
monic (blue dashed line) are also displayed for comparison. 

modulation and will naturally be present in the model (hence 
the smoothing); the second is due to the fact that appearing 
spots do not follow a smooth curve and present some disper¬ 
sion originating in the dynamo process. To take this contribu¬ 
tion into account, we calculated a second order polynomial 
dispersion d to the smoothed curve Rsmth. expressed as: 

d — d{j + ^ilRsmth "b ^2f^smth 

where da, d\, d 2 are derived empirically from observations. 
We then randomly added this dispersion d to the smoothed 
data Rsmth to obtain a daily input Wolf number R (Fig. E] 
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red dots, middle panel). The Wolf number (R, R) is a com¬ 
bination of the number of individual spots and of the number 
of spot groups. As the input of our simulation is the total 
number of spots N (which accounts both for the individual 
ones and for those in groups), we have derived a calibration 
to deduce this total number of spots from an activity level 
expressed in Wolf number units: 

N - Hi) + n\R (2) 

where hq, n\ are derived empirically from observations. 
Using Eq. |2] we obtained the theoretical total number of 
sunspots expected on the solar surface at each time step. The 
conversion no.i and dispersion £/o,i ,2 coefficients used in our 
simulations are given in Table. [1] By subtracting the number 
of sunspots already existing at this step, one can finally 
determine the number of new sunspots to be generated at 
this step in the simulation. 

We finally display in Fig. [T](bottom panel) the power spectra 
of the observed R and simulated input /? Wolf number for 
comparison purpose. The periodograms are quite different, 
especially at the rotational period and its first harmonic. This 
is expected as we do not introduce the rotational modulation 
(which is naturally present in the observed Wolf number) 
in our input Wolf number R but rather at the next step of 
the simulations (see above). Concerning the high-amplitude 
(and short-period) variations of the observed Wolf number 
in the high activity solar period, we consider that these 
variations are not due to a lack of spots at some moments 
in our model compared to the observations, but rather to 
the presence at these times of a group of large spots or an 
active cluster in the solar observations (see below for the 
spot size and spot filling factor distributions). Our model 
would require a more complex spot distribution function 
to reproduce such time-to-time appearances of large spot 
groups and we decided not to complicate it at this stage. 

We assumed an activity cycle length of 12.5 years (length of 
solar Cycle 23) and a time step of one day, leading to a total 
number of 4566 iterations or days. Furthermore, we added 
a random scattering of the time steps, with an amplitude 
of 4 hours around the regular daily period, to mimic real 
observations. 

2. Spatial and temporal distributions of the structures. At 
the beginning of the activity cycle, sunspots appear at 
medium latitudes, with some dispersion. As the activity 
level increases, the mean appearance latitude of sunspot 
slowly decreases in absolute value. This migration of active 
regions towards the equator continues during the whole 
cycle leading to a low mean appearance latitude at the end of 
the cycle. In our simulations, the mean latitudes of sunspot 
appearance at the beginning and at the end of the cycle 
were fixed at +22° and + 9°, respectively. These values are 
derived from the spot catalogs we used in Papers I and II. A 
linear law was then used to fit the decrease of the sunspot 
appearance latitude during the cycle, and we added a latitude 
scattering to reproduce the butterfly diagram. 

We also had the possibility to create an asymmetry be¬ 
tween the activity patterns in the northern and southern hemi¬ 
spheres, i.e. by adding more spots in one of them. As such a 
significant asymmetry has not been detected on the Sun dur¬ 


ing past solar cycles, we decided to put the same proportion 
of spots on each hemisphere. 

Spot groups are finally well-known not to appear at random 
longitudes, but around the so-called active longitudes. Such 
preferred longitudes h ave been detected on the Sun as well as 
on several other stos iBerdvugina & UsoskinL 120031 llvanovl 

l2007ilLanza et al.Ll20()9ll2010t) . Their origin is still partly un¬ 

explained and is probably due to the dynamo process. Some 
trends seem to emerge from observations : there are gener¬ 
ally two (seldom three) persistent active longitudes per hemi¬ 
sphere, shifted by 180°. They form a rigid structure, although 
they cannot be fix ed in a reference frame becaus e of the dif¬ 
ferential rotation (IBerdvugina & UsoskinL 1200^ . 

On the Sun, magnetic field and new activity seem to emerge 
at locations whe r e acti vity is already present (see e.g. 
iHarvev & Zwaanl Il993l) . Thus, instead of putting active 
longitudes with fixed longitude values in our model, we 
decided to build active longitudes from the position of 
already present sunspots. A fraction of new sunspots then 
appears in a restricted area in longitude around already 
present sunspots, while the remnant fraction is uniformly 
distributed in longitude. The fraction of new sunspots 
appearing in theses active areas and the longitude extension 
of the latter are reported in Table. [1] 

3. Large scale dynamics. We describe here the global motion of 
magnetic features at the solar surface during their lifetime. 
This motion can be decomposed in two components, namely 
a longitudinal and a latitudinal motion. 

The Sun is well known to exhibit a differential rotation in lat¬ 
itude: active structures rotate faster at the equator than at the 
poles. The longitudinal motion or rotation rate of the active 
structures a> (in degree per day) is therefore a function of the 
latitu de 6, according to the following equation (e.g. IWardL 
[T^lMeuni^ 12001 : 

CO = too + co\.sin^(0) + tO2.sin‘^(0) (3) 

The latitudinal motion or meridiona l flow M of magnetic fea¬ 
tures is derived from iKomm et al.l (Il993h . It is a poleward 
motion in each hemisphere, which is also a function of the 
latitude : 

M(6) = a.sin(28) + I3.sin(48) (4) 

4. Spot properties. Spots appear preferentially with a given 
size, increase in area rapidly before slowly decreasing. 
The spot growing phase is very quick (i n averag e 10 t o 
11 times faster than the decay phase, see iHowar Sum. 
It is then often shorter than our one-day timestep for the 
small and mid-sized sunspots and of a few days for the 
largest ones. As we focus here on long-period planets, and 
as the active structure growing phase is smaller or of the 
order of our simulation time step, we decided not to include 
it in our model for the moment and assumed that spots 
appear with their maximal size and then decrease during 
all their lifetime. We will however include a description of 
the growing phase in future works to better describe the 
variations at small timescales. Thus the spot distribution 
can be described with only the initial size distribution and 
a d ecay law. According to the La Laguna classification 
(see iMartinez Fillet et ^ 119931 and ref. therein), one can 
distinguish between isolated spots (La Laguna type 3) 
and spots belonging to complex groups (La Laguna type 
2). In compliance with this classification, we built two 
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Table 1. Input parameters. Free parameters are displayed in italic. 



Parameter 

Value 

Unit 

Reference 

Solar activity 

Reference cycle 

Solar Cycle 23 




Cycle length 

12.5 

[year] 



Time step 

1 

[day] 



Time step random dispersion 

4 

[hour] 



Wolf number normalization 

no = L32 





ni =0.148 




Wolf number random dispersion 

do = 6.922 





= 0.7594 





d2 = -0.00348 



Spatio-temporal 

Mean start latitude 

±22 

r] 


distribution 

Mean end latitude 

±9 

[°] 



Standard lat. dispersion 

6 

[“] 



Max. lat. dispersion 

20 

r] 



North/South asymmetry 

0.5 

- 



Active longitude spot fraction 

0.4 

- 



Active longitude extension area 

±20 

[“] 


Large scale dynamics 

Spot differential rotation 

cjo = 14.523 

[7day] 

Ward (1966) 



coi = —2.688 

[7day] 

This paper. 



Cl>2 — 0 

[7day] 

This paper. 


Facula and network 

cOfjQ = 14.562 

[7day] 

Meunier (2005) 


differential rotation 

= —2.04 

[7day] 

This paper. 



cJi,2 = -1.49 

[7day] 

This paper. 


Meridional flow, all structures 

a = 12.9 

[m s“*] 

Komm et al. (1993) 



1.4 

[ms“*] 

This paper. 


Stellar radius 

IR© = 696400 

[km] 


Spots properties 






Isolated spots 





Total fraction 

0.4 

- 

Martinez Fillet et al. (1993) 


Mean initial size 

46.51“ 

[pHem] 

Baumann & Solanki (2005)* 


Standard size deviation 

2.14 

[pHem] 

This paper. 


Max. size 

1500 

[pHem] 

Papers I and II 


Mean decay 

-18.9 

[pHem/day] 

Martinez Fillet et al. (1993)“ 


Median decay 

-14.8 

[pHem/day] 

This paper. 

Complex spot groups 


Total fraction 

0.6 

- 

Martinez Fillet et al. (1993) 


Mean initial size 

90.24“ 

[pHem] 

Baumann & Solanki (2005)“* 


Standard size deviation 

2.49 

[pHem] 

This paper. 


Max. size 

5000 

[pHem] 

Papers I and II 


Mean decay 

-41.3 

[pHem/day] 

Martinez Fillet et al. (1993)“ 


Median decay 

-30.9 

[pHem/day] 

This paper. 

Both spot types 


Min. decay value 

-3 

[pHem/day] 



Max. decay value 

-200 

[pHem/day] 



Min. spot size 

10 

[pHem] 

Papers I and II 

Faculae properties 

q (facula-to-spot ratio) 





Mean log(q) 

0.8 

- 



Standard deviation (log(q)) 

0.4 

- 



Min.-Max. log(q) 

0.1-5 

- 



Mean decay 

-27 

[pHem/day] 



Median decay 

-20 

[pHem/day] 



Min. facula size 

3 

[pHem] 

Papers I and II 

Network properties 

Diffusion coefficient 

300 

[km^s-‘] 

Schrijver (2001) 


Remainder fraction for decay 

0.975 

[/day] 



Min. size 

3 

[pHem] 

Papers I and II 


Facula fraction recovered in network 

0.8 

- 



Notes. For these two parameters, the value taken from iBaumann & Solankil (l200.'5h was multiplied by a factor 1.54 to better fit the spot catalogs 
used in Papers I and II. see the Table 1. of the paper, total area of single spots, snapshot method, Table 3., standard cycle. Table 1., total 
area, snapshot method. Table 1., standard cycle. 


distinct spot distributions. Both obey to similar evolution 
laws, but with different parameters. The initial spot size 


di stribution can be fitte d with a log-normal law, according 
to iBaumann & Solaris (l2005h . The two input parameters 
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are therefore the mean spot size and the standard spot size 
deviation. We took the parameters use d in our simulations 
from t he snapshot model developed in iBaumann & SolariO 
(l2005h. The spot decay l aw has been parametrized by 
iMartinez Fillet et aH (1 19931) as a log-normal law, with two 
input parameters ; the mean and median decay values. We 
also put an upper size limit for the two types of spots, and 
a lower threshold of 10 /rHem for all spots, chosen to be in 
good agreement with the observed spots distribution used in 
Paper II. 

5. Facula properties. On the Sun, spots and faculae form ac¬ 
tive regions, spots being surrounded by large faculae. The 
ratio between the surface covered by f aculae and spots (here¬ 
after the siz e ratio) has been studied (IChanman et all 1 19971 
l2Q0lLl2()Tll) . The authors found an average size ratio ranging 
between about 13 and 45 over solar cycles 22 and 23, de¬ 
pending on the method used to estimate the facular contrast. 
However, this corresponds to an average ratio measured for 
structures that can be at any state of evolution and any size. 
As for our simulation, we need instead an initial size ratio. 
Indeed, each time we add a spot, we add a facula at the same 
place. The facula size distribution then depends only on the 
initial size ratio. As there is little in the literature about this 
quantity, it is one of our few free parameters. We assume a 
log-normal law to describe the initial facula-to-spot size ra¬ 
tio. Facula decay is then described by two processes: 

- First the “classical” decay is described as for dark spots, 
with a log-normal law. This decay corresponds to a cer¬ 
tain surface being lost by the facula at each time step of 
the simulation. 

- Second, a given proportion of this facula lost surface is 
converted into network (the rest being completely lost, 
for example by flux cancellation or submergence): fac¬ 
ula fragments break away from active regions and diffuse 
over the surface. 

As the facula distribution is based on the spot distribution, 
the facula growing phase is not reproduced for the moment 
in our model. The facula growth and decay are longer 
than the spot on es in average, due to their longer lifetimes 
(I Howard Il992h . but th e facula growth is still much shorter 
than its decay (see e. s. iHoward 199 iT) . of a few days then. 
We decided not to reproduce it to remain coherent with 
the spot distribution and as we focus mainly on timescales 
between a fraction of the rotational period and a complete 
activity cycle, and on long-period planets. 

6. Network properties. We chose to make the network origi¬ 
nates in the decay of faculae (see point 5). Network behav¬ 
ior is then managed by its diffusion coefficient and by its 
decay rate. The network diffusion coefficient has been stud- 
ied in many papers w hich provide a great range of values 
(ICadavid et ai.Lll9^ and ref.therein). 

As these authors pointed out, there is a large gap between 
the diffusion coefficient values found by modeling (which 
are usually higher than 500 km^s ') and the values found 
with magnetograms analysis (which span from a few tens of 
km^ s“' to less than 300 km^ s“'). They also raised out the 
possibility for the diffusion coefficient to be time-dependent. 

W e finally d e cided to keep the value of 300 km^ s * used 
by ISchriiveij (1200 Ih as it stands well between the two ex¬ 
treme group s of values. We assum e the diffusion process to 
be isotropic (ICadavid et al.LlT99^ . As for the decay rate, we 
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Fig. 2. Upper panel, observed (black line) and simulated (red 
line) spot sizes. Lower panel: same for bright feature sizes. 

simply described it by fixing the network remainder fraction 
at each time step (it is one of the model free parameters). We 
finally put a minimal size of 3 /rHem for all bright structures 
(faculae or network), again in compliance with Paper II. 

2.3. Comparison with the observed activity pattern 

Here we compare the structures generated with our parameter¬ 
ized model to the observed structures used in Paper II. To do this, 
we use the spot and bright feature lists over Cycle 23 that we re¬ 
trieved from sunspot catalogs and MDI/SOHO magnetograms, 
respectively (Paper II). These data sets have some gaps, leading 
to a coverage of 3586 days, for a total duration of 4171 days, 
whereas our simulations have a length of 12.5 years with a daily 
sampling, i.e. 4566 days (length of Cycle 23). We therefore ap¬ 
plied to our outputs the same calendar as in Paper II so as to 
make a relevant comparison. The full time series will be studied 
in Sect.|4] 

2.3.1. Size distributions 

The size distributions (in /rHem) of our simulated structures are 
displayed in Fig. and compared with the observed ones from 
Paper II. We want to highlight the fact that very few of the input 
parameters detailed in Table [T] are free parameters. For most 
of them, we retrieved the input laws and parameters from the 
literature so as to build a fully parameterized model of solar-like 
activity. 

For the dark spot properties, we have no free parameters 
(i.e., they all originate from the literature). Remarkably, we find 
a good agreement between the spot size distribution simulated 
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Fig. 3. Upper panel: observed (left) and simulated (middle) spot latitude distribution over one Solar cycle. Right: Histogram of the observed 
(black) and simulated (red) spot latitude distributions. Lower panel: same for all bright features (faculae and network). Note that in the case of 
bright features, the histogram has its .t-axis in log-scale. 


this way and the observed one. The lack of a few large spots 
(about 40 in all) for the simulated distribution is explained by 
the fact that the log-normal law usually used in the litterature 
slightly underesti mates the number of very lar ge spots (see also 
Figs. 2 and 3 in I Baumann & Solan5 l2005h . The short-term 
differences reported above in the Wolf number distributions are 
then likely to originate in the time-to-time appearance of these 
very large spots (and not to a global lack of spots in our model, 
as the simulated distribution fits well the observed one). As for 
smaller spots, the observed distribution is discretized, which 
partly explained the gap (the remaining difference is again due 
to the fact that the observed distribution slightly differs from the 
log-normal law). This small discrepancy in the distribution for 
small spots will not influence significantly the RV signal as the 
RV effect of small spots is negligible compared to the larger 
ones. 

For the bright features (i.e., faculae and network, hereafter), 
there are three free parameters (the facula-to-spot initial size ra¬ 
tio distribution, the facula fraction recovered in network and the 
network decay rate). These parameters are poorly documented, 
so we decided to adjust them to better fit the observed distribu¬ 
tion. Overall, the simulated size distribution is in good agree¬ 
ment with the observed one (Fig. |2]). However, there are some 
remaining differences which are inherent to the model parame¬ 
ters. The size distribution for all observed bright features does 
not follow a log-normal law. For the simulated bright features, 
we remind that we chose to make the facula initial size distribu¬ 
tion directly dependent on the spot distribution, and the network 
originate in the decrease of the faculae. This explains thus the 
“flattened-S” shape of the simulated distribution, which corre¬ 
sponds to the addition of the facula and network respective dis¬ 
tributions. Reproducing even more accurately the shape of the 


observed distribution would require to make our model much 
more complex, for example by making the bright feature decay 
size-dependent or cycle-dependent. Given our goals, such a level 
of complexity is not necessary. As larger faculae have the largest 
influence on the RV, we decided to fit the distribution in priority 
for the highest sizes (i.e., for sizes greater than 5000 yuHem) by 
adjusting our free parameters. There is then still a discrepancy 
between the observed and simulated distributions for the small¬ 
est sizes, but we consider that its effect on our observables will 
be mostly negligible at our level of precision. 

2.3.2. Latitude distributions 

Another way to compare our simulated pattern to the observed 
one is to study the active structure latitude distributions. The lat¬ 
ter are displayed in Fig. [3] By doing so we will be able to validate 
the large scale behavior of our model. The latitude distributions 
of observed and simulated structures are indeed in good agree¬ 
ment. The kind of oscillations that we can distinguish in the ob¬ 
served bright feature latitude distribution at high latitudes and 
with a 1-year period comes from a “seasonal” effect. 

The solar observations are not made exactly in the plane of 
the solar equator, and the active structure distribution seems to 
be slightly shifted towards one visible solar hemisphere or the 
other, depending on the observation time, due to the noise level 
in the observed magnetogramfl This effect is also visible on 
the bright feature latitude histogram, where the observed distri¬ 
bution is spread a little more towards higher latitudes. 


^ As the noise level in the MDI magnetograms increases from the 
center to the edges of the images, the detectability level of a given struc¬ 
ture is different between the two hemispheres, depending on its latitude. 
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Projected spot filling factor 



Projected focuio filling foctor 


Fig. 4. Projected spot (top) and bright feature (bottom) filling factors (in yuHem), for observed (black) and simulated (red) structures. Left, projected 
filling factor time series (dots) and averaged over 30 days (solid lines). Note that a vertical offset has been put on the averaged curves for visibility. 
Middle: Rms of the projected filling factor over 30-day intervals. Right: histogram of the projected filling factor. 


2.3.3. Filling factors 

Moreover, we compare the filling factors of the projected active 
structures over the full solar cycle for the observed and simu¬ 
lated patterns in Fig. ID We find the observed and simulated fill¬ 
ing factor time series to be in very good agreement. The main 
difference between the observed and simulated distributions is 
located during the high-activity period; the activity peak present 
in the observations between approximately JD 2452200 and JD 
2452400 is not well reproduced in the simulations. This differ¬ 
ence does not come from a global lack of active structures in the 
high activity period of our model because i), the simulated daily 
spot number comes from Wolf number observations and ii), the 
facula distribution directly depends on the spot one (Sect. 12.21) . 
It rather originates in an occasional concentration of a few great 
structures in the observed pattern that is not reproduced in our 
model. Thus, this discrepancy is not a bias of our model but can 
be rather considered more as statistical noise. In the case of the 
bright features only, we also note a slight discrepancy between 
the observed and the simulated filling factors during the low ac¬ 
tivity period (the simulated filling factor being higher than the 
observed one). We attribute this to the larger number of very 
small bright features injected in our model (see before). 

Finally, we compare the ratio of the bright feature and spot 
filling factors over the solar cycle in Fig.|5]to estimate the evo¬ 
lution of the facula-to-spot size ratio. We find the ratios for the 
observed and simulated time series to be in very good agreement. 
One can also note that the size ratio reaches high values (from 50 
to 150) during the low activity period (i.e. at the beginning and 
at the end of the solar cycle), whereas it is much lower (around 
20) during the major part of the cycle. This is because during pe¬ 
riods of low activity, the number of dark spots is very low (equal 
or close to zero), whereas there is always a “background noise” 
due to bright features, and especially the network. 
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Fig. 5. Facula-to-spot size ratio averaged over 30 days, for ob¬ 
served (black) and simulated (red) structures. Blue solid line: 
average value over the cycle. 


We find the average value of the size ratios over the total cy¬ 
cle to be very close to each other, with values of 25.7 and 26.0 for 
the observed and simulated structur es, respectively. These values 
are in agreement with the studies of lChanman et al.l (Il997ll200ll 
201 ih. We also r emark that the size ratio time series displayed by 
Chapman et alJ (1201 Ih look very much like ours over cycle 23, 
with peaks above 100 during low activity periods and minima 
around 25 during high activity ones. 

Overall, we can conclude that our simulated solar activity 
pattern matches well the observed one and is thus a reliable 
model of a solar-like star magnetic activity. However, to defi¬ 
nitely establish its validity, we compare the resulting RV and 
photometric time series. In particular, we check the impact of 
the differences between the two patterns on the time series. 
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Fig. 6. Left, observed TSI from iFrohlich & Lean! (Il998l) (black) and reconstructed photometric contribution from all simulated struc¬ 
tures (red). The photometric time series (dots) are displayed as well as their average over 30 days (solid line), with a vertical offset 
on the latter for visibility. Middle: Rms of the photometric time series over 30-day intervals. Right: Histograms of the photometric 
time series. 


3. Comparison with the Sun observations 

In this section, we first describe the way we simulate the time 
series for the different observables and how we take into account 
the stellar properties (limb-darkening, active structures tempera¬ 
ture contrasts, impact on convective blueshift, physical and geo¬ 
metrical stellar properties). We then compare the RV time series 
simulated in the case of an edge-on solar-like star to the RV time 
series obtained with the observed activity pattern described in 
Paper II to establish the validity of our model. 

3.1. Description of the simuiations 

3.1.1. Limb-darkening 

In our previous papers, the center-to-limb darkening was linear 
with respect to p (with/r - cos(6), ff being the angle to the center 
of the solar disc. lDesort et al.Ll2007h . as follows: 

I(p)^l-e + ep (5) 

with e = 0.6. In this case, the limb-darkening was not 
temperature-dependent and therefore applied to both the inac¬ 
tive photosphere and the active structures indiscriminately. Here, 
we decided to change our limb-darkening law, so as to: i) have a 
more accurate one; and ii) use a temperature-dependent one, that 
could be adaptable to different types of stars. We took the non¬ 
linear limb-darkening law from IClaret & Hauschildd (l2003h : 

§=i-|:«(i-.h (6) 

where 7(1) stands for the intensity at the center of the stel¬ 
lar disc. The four Claret limb-darkening coefficients Ok are 
the bolometric coefficient s taken from ATLAS models (see 
IClaret & Hauschild3.l20Q.3h . These coefficients are temperature- 
dependent. That is why such a law is better adapted to our simu¬ 
lation tool, as the effective temperature is one of the input param¬ 
eters, and such a law can be extrapolated not only for the Sun, but 
for any star for which the effective temperature is known. It also 
means that the limb-darkening coefficients applied to the spot¬ 
ted stellar surface will be slightly different from the ones applied 
to the inactive photosphere. We therefore use a different limb- 
darkening for the photosphere and for the dark spots, depending 
on their respective temperature. As faculae show a very strong 
contrast variation depending on their position on the disc due to 
more complex processes, we directly define their contrast Cpi(p) 


with respect to the limb-darkened photosphere (see below) with¬ 
out using the Claret law according to their temperature. 

3.1.2. Spot and bright feature contrast and photometric time 
series 

To estimate the contrast of an active structure (dark spot or bright 
feature) with respect to the quiet photosphere, we use the proce¬ 
dure described in Paper II. This procedure is an independent, 
preliminary step to our simulations (i.e., to the generation of the 
spectra and of the RV and photometric time series), and is aimed 
at determining the active structure contrasts that we will use as 
input parameters. Starting from a given range of values for the 
spot and bright features contrasts, and based on observed struc¬ 
ture patterns, we build (for each set of contrast values) two time 
series representing the respective relative contributions of dark 
and bright active structures to the stellar irradiance. A min¬ 
imization between the sum of these contributions (quiet photo¬ 
sphere, dark and bright features) and the observed total solar ir¬ 
radiance (hereafter TSI) taken from iFrohlich & Lead (Il998l) is 
then performed over the solar cycle 23. We adopt the same value 
for the quiet Sun reference as in Paper II, i.e. 13 65.46 W.trT^, 
which is very close to the average value found by ICrouch et al.l 
(l2008h with their TSI model over twelve solar cycles. Contrary 
to Paper II, we include in the procedure the influence of the 
center-to-limb darkening, so as to take into account the Claret 
limb-darkening law we will now use in our simulations. We end 
up with a spot temperature deficit ATsp - -6Q5K and a facula 
contrast Cp, = 0.131618 - 0.218744yu H- 0.104757//2. We then 
use these contrasts as our input parameters in the simulations. 
As in Paper II, we finally compare the sum of the photometric 
contributions of spots and bright features obtained with ou r sim¬ 
ulations with the observed TSI of iFrohlich & LeanI (Il998h . so as 
to check the validity of our activity model. The comparison is 
done on 2263 points, ranging from 1996 to 2003. 

The irradi ance obtained wi t h our simulations and the ob¬ 
served TSI of iFrohlich & LeanI (Il998l) are displayed in Fig. |6] 
We match quite well the observed irradiance, excepted for some 
peaks in the high activity period that are not reproduced in the 
model. This can be explained by the fact that our model does 
not reproduce well the occasional appearance of very large ac- 

^ As in Paper II, the spot contrast is a temperature contrast, while the 
bright structure (facula or network) contrast is defined as C„; = -4——, 
with S pi and Sp/, the facula/network and quiet photospheric irradiances, 
respectively. 
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tive structures or large active structure clusters (see Sects. ISjTI 
and I2.3.3 T i. Consequently, the distribution of the simulated TSI 
is slightly narrower than the observed TSI one, and the averaged 
rms of the simulated TSI is about 15% lower in comparison (see 
Table O. However, the temporal evolution of the simulated TSI 
rms matches well the observed one. As already stated, we con¬ 
sider that these differences do not affect the signihcancy of our 
simulations as a model. 

3.1.3. Attenuation of the convective blueshift 

For all following simulations, we adopt the same value for the at¬ 
tenuation of the convective blueshift as in Paper II, i.e. 190 m s“*. 
We refer to this paper for justification. 

3.1.4. Building the spectra 

W e build the spectra from the input structure lists as described 
in iDesort et al.l (l2007l) and in Papers I and II. We assume a stel¬ 
lar mass of 1 Mq, a temperature of 5800 K and a rotational 
velocity of 1.9 km s ' at the equator. We keep in this section the 
stellar inclination i to 90°, corresponding to a star seen edge- 
on. Brie f ly, we use a synthetic spectrum from Kurucz models 
(lKuruczIll993h corresponding to a G2V star and apply it to each 
cell of the visible stellar 3D hemisphere divided into a grid. The 
spectrum is shifted to the cell radial velocity. In case of the pres¬ 
ence of an active structure, it is weighted with a black-body law, 
taking into account the active structure contrast with respect to 
the quiet photosphere. To include the effect of the attenuation of 
the convective blueshift, the cell spectrum is redshifted by 190 
ms *. T his attenuation has oth er effects on the spectral lines, 
however bumusaue et al.l (1201 4l) show that only considering the 
RV shift is sufficient to estimate the RV effect of active regions. 

All cell spectra are then balanced by the cell’s projected sur¬ 
face and limb-darkening. We hnally sum up all cell contributions 
to obtain the stellar spectrum. 

3.1.5. Computation of the RV time series 

As in Paper II, we compute the RV using our Software for the 
An alysis of the Fourier I nterspectrum Radial velocities (SAFIR, 
see iGalland et al.l l2005h on the built spectra as if they were ac¬ 
tual observed spectra. We use only the wavelength range corre¬ 
sponding to the order jj 31 of the HARPS spectrograph, as done 
in Papers I and II. We obtain three distinct time series: two due 
to the respective photometric contributions of dark spots and 
bright features (hereafter the spot and facula time series, respec¬ 
tively); and the third due to the partial inhibition of the convec¬ 
tive blueshift in the active structures (hereafter the convection 
time serie). For the latter, we consider only the inhibition of the 
convective blueshift (i.e. this effect is not weighted by the active 
region flux). We sum the three time series to obtain the total RV 
variations. Here we consider this sum to be a good approxima¬ 
tion of the real RV variations given that the convective blueshift 
is dominant (see below). 

3.2. Validation of the activity model 

3.2.1. RV time series basecd on the observed patterns 

In this section, we first check if the change of limb-darkening 
law and the corresponding change of the active structure con¬ 
trasts have a signihcant impact on the resulting RV time series. 
To do so, we compute new RV time series corresponding to the 
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Fig. 7. RV time series computed with Claret limb-darkening and 
corresponding structure contrasts versus RV time series com¬ 
puted with linear limb-darkening (as in Paper II). Top: contri¬ 
bution of spots and faculae. Bottom: sum of all contributions. 


observed solar activity patterns used in Paper II, but this time 
with the new limb-darkening law and structure contrasts. We 
then compare these new RV time series to the ones obtained in 
Paper II, using the same temporal sampling. We hnd the RV time 
series to be closely correlated and thus in good agreement. We 
display in Fig. [T] the RV time series computed with the Claret 
limb-darkening versus the RV time series computed with the lin¬ 
ear limb-darkening. We hnd both the slope of the hts and the 
correlation of the RV time series to be close to 1, as well in the 
case of the spot-i-facula RV time series as for the total RV ones. 
Despite the simpler limb-darkening law and a different active 
structure contrast, the results shown in Paper II are still valid for 
estimating the RV effect. 

3.3. Comparison between the RV time series based on the 
observed and simulated patterns 

We compare here the RV time series based on the observed and 
simulated activity patterns (both being computed with the Claret 
limb-darkening) so as to assess the validity of our activity model. 
We compare the RV time series for the different contributions 
(spot, facula, convection and all). All time series, as well as their 
dispersion over the cycle and their histograms, are displayed in 
Fig. [8] We hnd the amplitudes of the observed and simulated RV 
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Fig. 8. From top to bottom: Spot, facula, convection and total contribution to the RV time series, reconstructed from observations 
{black] same activity pattern as in Paper II) and simulated with our model (red). Left: RV time series {dots). The timescale is the one 
from Paper II. For the convection and total contributions, the RV time series averaged over 30 days are also displayed {solid line), 
with an offset for better visibility. Middle: Rms of the RV over 30-day intervals. Right: histograms of the RVs. 


time series to be in very good agreement for the spot and fac¬ 
ula time series, with very similar histograms. As for the RV time 
series corresponding to the convective component, it is closely 
related to the facula filling factor (with a Pearson correlation co¬ 
efficient of 0.97) and widely dominates the total RV signal, as 
in Paper II. The main visible difference comes from the activ¬ 
ity peak at around JD 2452200-2452400 in the RV reconstructed 
from observations, which is not echoed in the simulated RV. We 
already discussed the origin of this difference in Sect. 12.31 In 
the low activity period, the averaged amplitude of the simulated 
convection time series is about 20% higher than for the observed 
one (see Table[3]|. We attribute this small discrepancy to the ex¬ 


cess of very small bright features in the model that we discuss 
Sectl23] 

To have a better understanding of the RV signature of our 
simulated activity pattern, we also display the RV rms computed 
over 30-day intervals over the cycle in Fig. 0 This gives an idea 
of the temporal evolution of the RV time series dispersion. We 
hnd the evolution of the RV dispersion during the solar cycle to 
be in good agreement for the observed and simulated time series. 
We also provide the RV rms for the different components taken 
over the complete cycle and for low and high activity periods in 
Table |2l and the RV amplitudes in Table |3 The low and high 
activity periods are the same as in Paper II. We note that: 
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Fig. 9. Lomb-Scargle periodograms of the RV time series based 
respectively on the observed (black solid line) and simulated 
(red solid line) activity patterns, alongside with the 1% false- 
alarm probabilities (FAP, dashed lines). Top: spot-i-facula contri¬ 
bution. Bottom: sum of all contributions. Note that in both pan¬ 
els the periodograms corresponding to the simulations are verti¬ 
cally shifted for visibility. The solar rotation period and its two 
first harmonics are displayed in blue (solid, dashed and dashed- 
dotted blue lines). 

- For the spot and facula simulated time series, the rms is 
slightly higher than for the observed time series in the low 
activity period (about two times higher in the case of bright 
features). However it is not the case for the convective time 
series. 

- In the case of the convective component, the rms over the 
complete cycle is about 24% lower for the simulated data 
than for the observed data. When looking at the temporal 
evolution of the RV rms displayed Fig.|8] we can see that this 
discrepancy mostly comes from a few high activity peaks in 
the observed time series. 

We finally compare the Lomb-Scargle periodograms of the RV 
time series based on the observed and simulated activity pat¬ 
terns. The periodograms are displayed in Fig. 0in the case of 
the “photometric” contributions (spots and faculae) and of the 
sum of all contributions. The respective periodograms for the 
observed and simulated patterns remarkably present the same 
characteristics: 

- For the “photometric” component, the power is mostly con¬ 
centrated at the solar rotation period and its two first harmon¬ 
ics. 


Table 2. RV rms (in m/s) for the different components. 


Simulated pattern 

spots 

faculae 

sp+fac 

conv. 

total 

All 

0.34 

0.32 

0.33 

1.98 

2.00 

Low 

0.17 

0.16 

0.14 

0.48 

0.51 

High 

0.47 

0.38 

0.47 

1.36 

1.52 

Observed pattern^ 

spots 

faculae 

sp-l-fac 

conv. 

total 

All 

0.37 

0.25 

0.32 

2.59 

2.62 

Low 

0.10 

0.08 

0.07 

0.47 

0.47 

High 

0.52 

0.35 

0.42 

1.50 

1.53 


Notes. i.e. structure pattern from Paper II and Claret non-linear limb- 
darkening. 


Table 3. RV rms and peak-to-peak amplitude (in m/s) and relative rms 
photometry for the three periods. 


Simulated pattern 

rms RV 

ampl RV 

rms phot 

All 

2.00 

11.3 

2.83 10^^ 

Low 

0.51 

2.5 

1.06 10-“* 

High 

1.52 

7.7 

2.97 10-“* 

Observed pattern^ 

rms RV 

ampl RV 

rms phot 

All 

2.62 

11.4 

3.6 10-“ 

Low 

0.47 

2.1 

1.2 10-“ 

High 

1.53 

8.4 

4.5 10-“ 


Notes. i.e. structure pattern from Paper II and Claret non-linear limb- 
darkening. 


- For the sum of all contributions, there is still power at the ro¬ 
tation period, but there is much more power at longer periods 
(between 300 and 1000 days) due to the long-term activity 
cycle. 

Nonetheless, we note a significant difference: for the 
spotH-facula RV signal, the peaks at half the rotation period are 
emphasized compared to the peaks at the rotation period in the 
case of the observed pattern, whereas it is the contrary for our 
simulated pattern. We already noticed and discussed in Paper 
II the predominance of the power at half the rotation period in 
the RV time series derived from observed solar patterns. We ex¬ 
plained it by the presence of two symmetrically active longitudes 
on the solar surface. Since here we do not impose two active lon¬ 
gitudes separated by 180° over the whole time series but rather 
variable active longitudes, this effect may not be very important. 
An alternative explanation could be that the difference between 
the periodograms corresponding to observed and simulated pat¬ 
terns may originate in the active structure quick growing phase 
(which is at present not reproduced in our model). 

This may also explain the smaller RV and photometric rms 
when compared to the observations. 

As detailed above, small discrepancies in amplitude and dis¬ 
persion on short timescales are present both in the photometry 
and in the RV between the time series derived from observed 
and simulated activity patterns. We consider that these differ¬ 
ences likely originate in two sources: first, the spot size distribu¬ 
tion, where the few larger spots are hard to model; and second, 
the active structure growing phase. Overall, we consider that the 
temporal evolution over the complete cycle of both the amplitude 
and dispersion of the simulated time series match well the time 
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Fig. 10. Left, projected spot filling factor at maximum activity (in 
fraction of the stellar surface) versus stellar inclination i. Right: 
same for bright features. 

series derived from the observed activity pattern. Hence, we con¬ 
clude that: i) these short-term discrepancies are not significant 
enough to question our model reliability; and ii) the comparison 
between the RV time series derived from observed and simulated 
patterns then assess the overall validity of our model. 

4. Photometric and RV time series of inciined 
soiar-type stars 

For solar-type stars of the same age as the Sun, it is commonly 
assumed that active structures are mostly concentrated in a belt 
around the stellar equator (even if bright features are much more 
dispersed in latitude). A different inclination of the stellar ro¬ 
tation axis should then have a significant effect on the various 
time series that needs to be investigated. In this section, we per¬ 
form the same simulations as above, but for different inclina¬ 
tions of the stellar rotation axis. We consider inclinations be¬ 
tween i - 10° (i.e., for a star seen nearly pole-on) and i - 90° 
(star seen edge-on), with a sampling of 10°. The time series are 
studied with their original temporal sampling, i.e. 4566 days and 
no gaps. 

4. 1. Photometric time series 

The impact of the inclination i of the stellar rotation axis 
on the long- t erm s olar irradiance variations was studied by 
iKnaack et alJ (l200lll . The authors computed the solar irradiance 
corresponding to a 3-component model (quiet Sun, dark spot 
and bright facula) at activity extrema and and for a variable 
inclination. As for their active region distributions, they used 


Fig. 11. Top: Relative variation of the TSI over the activity cycle (i.e., 
peak-to-peak amplitude of the TSI time series averaged over 30-day 
intervals, expressed in fraction of the quiet Sun irradiance reference) 
versus inclination i (solid line). The respective contributions of spots 
(stars) and bright features (diamonds) are also displayed. Bottom: short¬ 
term TSI dispersion (rms of the residuals after substraction of a 30-day 
averaged time series) versus i. 


simple active latitude belts with contrast corresponding to the 
given active structure (i.e. no individual structure was intro¬ 
duced in their model). According to them, the apparent active 
structure surface coverage decreases with a decreasing i. They 
nonetheless expected an increase of the TSI with a decreasing 
i since bright features are limb-brightened while the contrast 
of dark spots is roughly independent of the limb-darkening. 
Thus, when going from an edge-on toward a pole-on configu¬ 
ration, the impact of dark spots would decrease following their 
apparent covered surface. On the contrary, the decrease of the 
apparent surface covered by bright features would be at least 
compensated by their increased contrast since we would see 
them mainly on the limb. We display in Fig. [10] the evolution 
of the spot and facula projected filling factors with i. We find 
the projected spot filling factor to decrease by about 35% when 
going from edge-on (i- 90°) to nearly pole-on (i= 10°), and the 
projected facula filling factor to decrease by about 43% . Thes e 
results are quite similar to those found bv iKnaack et al.l (l200lh . 
We note that the filling factor does not tend towards 0 for i ^0°. 

The relative variation of the TSI during the cycle (which is 
of the order of 0.1% of the quiet Sun irradiance) vs. i is dis¬ 
played in Fig, nn (upper panel). We find it to increase by only 
14% when going from i- 90° to /= 10°. This is much s maller 
than the 40 + 10 % increase predicted bv iKnaack et al.l (12001 h 
with comparable input parameters but a much simpler model. 
Yet, it points towards the same verdict as pulled through by 
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iKnaack et al.l (1200 ih . i.e. that photometric variations of seem¬ 
ingly inactive Sun-like stars cannot be explained by an inclina¬ 
tion effect. We conclude that for a solar-like star with a simi¬ 
lar activity configuration, the effect of an inclined rotation axis 
on the long-term variations of the total iri'adiance (of the order 
of the solar cycle length) is relatively small. On the contrary, 
the TSI short-term v ariations (which were indeed not studied by 
iKnaack et al.Ll200Th are strongly impacted by the inclination. As 
illustrated in Fig. [TT] (lower panel), we hnd the TSI short-term 
dispersion to be decreased by a factor ~ 6 when going from 
i- 90° to i- 10°. A possible explanation for such a decrease is 
the following: 

- First, for smaller inclinations, we mainly see the effect of 
bright features as they are more extended towards the higher 
latitudes than the dark spots. 

- Then, due to the activity configuration (where active struc¬ 
tures are mainly located on two belts on both sides of the so¬ 
lar equator), for small inclinations we see the same structures 
during all the rotation period and not during half a period. 

Therefore, in the case of a star seen nearly pole-on, the short¬ 
term dispersion of the irradiance (of the order of a few rotation 
periods) originates no more in the structure crossing of the vis¬ 
ible hemisphere for each half rotation period. It originates only 
in the structures appearance and decay. To confirm it, we display 
in Fig.[T2]the Lomb-Scargle periodograms of the simulated TSI 
for representative inclinations. For a star seen edge-on, the peri- 
odogram is dominated by power at the rotation period of the star. 
When going towards smaller inclinations, the signal at the ro¬ 
tation period gradually decreases until it disappears completely 
for a star seen nearly pole-on. On the contrary, the signal at a 
much longer period (which is likely induced by cycle-related 
long-term periodicities of the order of the cycle length) becomes 
increasingly preponderant for smaller inclinations. A reason for 
which the signals at long term periods (in the 600 to 1500-day 
range) become increasingly dominant with decreasing i in the 
TSI (and to a lesser degree in the total RV) periodograms could 
be the following: for nearly pole-on conhgurations, we see the 
active structures on one stellar hemisphere only, whereas we see 
them on two hemispheres for configurations closer to edge-on. 
This may induce a kind of an averaging effect on the long-term 
signal for the edge-on configuration and explain its increase with 
decreasing i. We hnally display in Fig.[T3]the simulated TSI for 
representative configurations. As we found above, the long-term 
amplitudes of the two time series are nearly the same; on the 
contrary the dispersion is widely reduced with i. 

4.2. RV time series 

We now study the impact of stellar inclination i on our simu¬ 
lated RV time series. Our first main result is that in contrast with 
the photometry, both the amplitude and the dispersion of the to¬ 
tal RV decrease with a decreasing i. This is well illustrated in 
Fig.[l4]where we display the RV time series for all activity com¬ 
ponents and for signihcant inclinations. The amplitude as well as 
the dispersion (RV rms) of the total RV taken over the complete 
cycle are decreased by a factor ~ 6 when going from an edge-on 
to a nearly pole-on conhguration. 

We characterize both the long-term and short-term variations 
of the RV signals to investigate deeper the impact of inclination. 
Our results are illustrated in Fig. [15] The “long-term” variations 
simply cori'espond to the signal taken over the complete stellar 
cycle. To study the short-term variations of the signal, we per¬ 
form a running average of the RV time series with a smoothing 


window of 30 days and substract it to the original data. Then the 
long-term variations should primarily be affected by the global 
activity cycle, whereas the short-term variations will come from 
rotation-related effects. 

1. Peak-to-peak amplitude: For the total RV and when going 
from i - 90° to i - 10°, the amplitude decreases by ~ 80% 
over the cycle and by ~ 85% on the short term. When tak¬ 
ing only the “photometric contribution” (i.e., spots and facu- 
lae) into account, the peak-to-peak amplitude is decreased by 
nearly 90% on the long-term, and by 87% on the short-term.. 

2. Dispersion: the total RV rms decreases by 70% with incli¬ 
nation on the long term, and by 85% on the short term. In 
the case of spots and faculae only, the decrease is the same 
on both timescales and is of nearly 85%. For a star with 
a solar-like activity pattern seen almost pole-on, a jitter of 
~ 0.5 ms“' can be expected if the observation timescale is 
of the order of the activity cycle length, and a jitter of ~ 0.2 
m s ' can be expected over a month. 

3. Ratio between “photometric” and convective components: 
we also study the evolution of the relative contribution of the 
“photometric” (i.e., due to spots and faculae) component to 
the total RV signal with the inclination. Over the cycle, the 
“photometric” fraction of the RV is reduced by a factor ^ 2 
when going from i - 90° to i = 10° (i.e., the convective 
component is increasingly preponderant with a decreasing 
inclination). On the contrary, on the short-term the “photo¬ 
metric” relative contribution remains nearly constant with /, 
at a level of ~ 0.3. This value is at least two times larger than 
the value of the “photometric” fraction on the long-term (i.e., 
~ 0.16 in the edge-on conhguration). This conhrms that the 
spot-Hfacula component has mostly a short-term effect on the 
RV signal, of the order of the rotation period. 

We hnally study the Lomb-Scargle periodograms of the RV 
time series for representative inclinations, hrst in the case of 
the sum of the hux contributions of dark and bright structures, 
and then when adding the convective component. The peri¬ 
odograms are displayed in Fig. [T2| Starting from i = 90°, the 
RV periodogram for the spot-i-facula signal is dominated by the 
stellar rotation. The power is concentrated at the rotation period 
and its two hrst harmonics, with about two times and six times 
less power for the hrst and second harmonics, respectively. 
The power is induced by the active structure crossing of the 
visible hemisphere during about half the rotation period (as 
the active structures are mainly located near and around the 
stellar equator). When i decreases, we begin to see the active 
structures on a time longer than half the rotation period, until 
we see them permanently for the nearly pole-on conhguration. 
There is thus nearly no remaining component in the spot-i-facula 
RV periodogram for i - 10°. This is in good agreement with the 
90% decrease in jitter already described. 

Starting from the edge-on conhguration, the total RV peri¬ 
odogram is dominated: i) by a serie of peaks at long period (200 
- 1000 days) due to the long-term activity cycle; ii) by peaks at 
the stellar rotation period and its three hrst harmonics. 

Expectedly, the peaks corresponding to the rotation pe¬ 
riod and its harmonics gradually decrease when going towards 
smaller inclinations, until they nearly disappear for i = 10°. 
We also note that the harmonics of the rotation signal decrease 
signihcantly faster than the rotation signal itself. On the con¬ 
trary, the long-term signal corresponding to the activity cycle 
decreases slowly when i decreases, but remains widely above 
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Fig. 12. Top: Periodogram of the TSI for significant inclinations (from left to right: i = 10°, 30°, 60° and 90°). The 1% false-alarm prohahilities (FAP) are displayed (blue dashed line), as well as 
the equatorial stellar rotation period and its two first harmonics (red solid, dashed, and dotted-dashed lines, respectively), and the 480-day period for the computing of the detection limits (green 
solid line). Middle: the same for the RV “photometric” component (spot-l-facula). Bottom: the same for the total RV. 
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Fig. 13. Simulated TSI for: i - 90° (black), i = 60° (purple), i = 30° 
(red) and i = 10° (orange). Vertical offsets have been introduced for 
visibility. 


Fig. 14. Simulated total RV for significant inclinations. From top 
to bottom: i — 90° (black), i - 60° (purple), i - 30° (red) and i - 
10° (orange). Vertical offsets have been introduced for visibility. 


the 1% false-alarm probabilities (FAP) for i = 10°. In the edge- 
on configuration, the power is equally distributed between the 
rotation and the long-term signals, but the latter become pre¬ 
ponderant for smaller inclinations. This is in agreement with the 
decrease of the ratio of the RV “photometric” component over 
the total RV signal on the long-term we show above. We con¬ 
clude that in the case of a star with a solar-type activity pat¬ 
tern, the stellar inclination strongly impacts both the observed 
stellar irradiance, RV variations and active regions filling fac¬ 
tor (and hence the observed chromospheric activity). Indeed, the 
short-term variations of these various observables strongly de¬ 
crease with a decreasing inclination. Then, on short to medium 
timescales (from a few days to a year, or well under the activity 
cycle length), it will probably be not possible to distinguish be¬ 
tween a solar-like star seen nearly pole-on and an inactive star. 
Provided that the time baseline is sufficient, on timescales of the 
order of the cycle length, the activity cycle should however still 
be detectable for small inclinations, but with a reduced ampli¬ 
tude. It would then not be possible to distinguish it from a nearly 
inactive star with a low-amplitude activity cycle. 

4.3. Potential impact of steiiar inclination on various 
correlations 

4.3.1. Correlation between total RV signal and facula filling 
factor 

In Paper IV, we used the solar Calcium (Ca) index (S-index or 
log(Rjjj.)) as a tool to correct partially the activity-induced RV 
from its convective component (using the observed solar activ¬ 
ity pattern). Indeed, the Ca index is correlated with the chro¬ 
mospheric plage fill ing factor and hence with the photospheric 
facula filling factor dShaniro et all l20l4 found a linear depen¬ 
dence between the Ca index and the facula filling factor in the 
case of the Sun). 


As already explained, in the case of a solar-like star seen 
edge-on, the RV signal is closely correlated to the facula fill¬ 
ing factor (as the convective component dominates the RV vari¬ 
ations, see Sect. l3.3l and Figs. When going toward smaller 
inclinations, we note that the correlation gets even stronger, with 
a Pearson correlation coefficient going from 0.97 for i - 90° to 
0.99 for i - 10°. This is in agreement with the fact that the rela¬ 
tive contribution of the spot-i-facula component to the RV signal 
decreases with i. It would finally mean that the Ca index correc¬ 
tion method can be used regardless of the stellar inclination. 


4.3.2. RV jitter, Ca Index and facula filling factor 


In recent studie s dWrighli 1200.51 ISantos et al. . 1201 Oi 


llsaacson & Fischeil l2010l iHillenbrand et al. l2014li . the 


average level of RV jitter was commonly taken as a proxy for 
the stellar magnetic activity level, along with the Ca index. 
For different samples of FGK stars, the authors found loose 
correlations between the average RV jitter and the mean Ca 
index, generally with a significant amount of dispersion. In 
our simulations, we find that the ratio of the mean facula 
filling factor over the RV rms (calculated over the complete 
activity cycle) is not constant over the range of inclinations i we 
explored. We display in Fig. [16] the facula filling factor to RV 
amplitude and to RV rms ratios. We find both ratios to increase 
towards smaller inclinations (in agreement with Sects. 14.11 
and 14.21 where the decrease of the RV amplitude and rms with i 
is more pronounced that the decrease of the facula filling factor 
with /). Provided that the Ca index evolves in the same way as 
the facula filling factor with i, this could partly explain the large 
dispersion in the (RV rms, mean Ca index) relation found in the 
studies cited above (as the stellar inclination remains unknown 
for most of the observed targets). It also means that a clear (RV 
rms, Ca) relation will remain hard to establish for future studies, 
unless the uncertainty on the stellar inclination can be removed. 
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Fig. 15. Characterization of the RV signal versus sin(i). Left. Long-term variations of the signal (RV time series as produced hy our simulations). 
Right. Short-term variations of the RV signal (residuals of the original RV data minus the 30-day averaged RV). From top to bottom: Peak-to-peak 
amplitude of the RV signal, RV dispersion (rms), ratio of the amplitudes of the spot-l-facula RV to the total RV, and ratio of the rms of these two 
signals as a function of sin(i). On the two top panels, the spot-l-facula {dotted line), convective (dashed line) and total RV (solid line) are also 
displayed as a function of sin(i). 
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Fig. 16. Top panel: Ratio of the mean facula filling factor over 
the total RV amplitude versus i. Bottom panel: Ratio of the mean 
facula filling factor over the total RV rms versus i. 


4.4. Detection limits 

In this section, we compute the detection limits for the total RV 
time series over the range of stellar inclinations explored above. 

4.4.1. Approach 

We use two different methods: the correlation-based method 
and the local power amplitude (hereafter LPA) method. Both 
were described in de tails and tested on real stellar RV data in 
iMeunier et al.l (l2012b . In brief, the first method makes the cor¬ 
relation between the periodograms of a generated planetary RV 
signal (with the same temporal sampling as of the data) and of 
the actual RV data to which the planetary RV signal has been 
added (i.e., it determines the correlation of the power of the (stel¬ 
lar data -H fake planet) periodogram vs. the power of the fake 
planet periodogram). The detection limit corresponds then to the 
minimal mass and period for which the correlation values are all 
above a given threshold, for 100 realizations spanning all orbital 
phases. The threshold (spanning here from 0.003 for i - 10° to 
0.05 for i - 90°) corresponds to the maximum of the correla¬ 
tions obtained for a very low mass planet (here ~0.6 MEauh- As 
for the LPA method, we compare the periodograms of the ac¬ 
tual RV data and of a given generated planetary signal (for the 
same temporal sampling), but within a localized period range 
around the given planetary period. The detection limit at this pe¬ 
riod corresponds then to the minimal mass for which the maxi¬ 
mum power of the planetary RV periodogram (within the limited 
period range) is always above the maximum power of the actual 
RV data periodogram in the period range. 

We compute here our detection limits for only one period of 
480.1 days, corresponding to a separation of about 1.2 au, as in 
Paper II. It is a representative value of the outer boundary of the 


HZ for a solar-type star. The total time span is always fixed to 
the complete simulation time span (i.e. 4566 days). As in Papers 
I and II, we assume the star to be observed eight months a year 
(meaning that we take into account only 2978 RV points instead 
of 4566). We assume different temporal samplings: all points 
(i.e. one-day sampling or 1:1), 1 point every 4 nights (1:4), 1 
point every 8 nights (1:8) or 1 point every 20 nights (1:20); that 
is to say that the (1:20) series has 20 times less points that the 
(1:1). We also assume different precisions on the RV (no added 
noise, 0.01, 0.05 and 0.1ms * noise levels), as we did in Paper 

l. 

Finally, we compute the detection limits vs. inclination for 
two different cases: 

1. In the first case, we compute the detection limits considering 
that the hypothetical planet is always seen orbiting edge-on 
for all stellar inclinations (i.e., there is an increasing spin- 
orbit misalignment with decreasing inclination). 

2. In the second case, we consider that the hypothetic planet al¬ 
ways orbits in the stellar equatorial plane (i.e. there is always 
spin-orbit alignment). 

Considering these two distinct cases is important as in most 
cases when searching for planets with RV, we do not know nei¬ 
ther the stellar inclination nor the inclination of the planet orbit 
with respect to the line of sight. As for the inclination of the 
planet orbit, it is taken into account by giving minimal masses 

m. sin(i) (with i denoting here the inclination of the planet orbit 
with respect to the line of sight) for detected planets. However, 
the detection limits are generally computed without knowing 
the stellar inclination itself, i.e. considering that the star is seen 
edge-on and that the activity-induced jitter is not reduced due to 
the projection effect. The detection limits we compute in case 1 
correspond to this configuration and thus to the “best-possible” 
detection limits. As we go toward smaller stellar inclinations, 
they also correspond to more and more unlikely orbital configu¬ 
rations. Indeed, most of the exoplanets found so far indeed orbit 
in or near the stellar equatorial plan^l (even if in the specific case 
of Hot Jupiters, a si gnificant fraction of systems show spin-o rbit 
misalignments, e.g. l Albrecht et al.Ll2012llTriaud et ani2014l) . 

On the contrary, the detection limits that we compute in case 
2, i.e. when considering a spin-orbit alignment, correspond to a 
more conservative but more realistic assumption. The detection 
limits for the correlation-based and LPA methods are displayed 
in Fig. [n] 

4.4.2. Results 

Comparison with previous results - For the edge-on config¬ 
uration, we can compare the detection limits (for the 1:1 sam¬ 
pling) to the detection limits computed in Paper i.e. on RV 
time series derived from observed solar activity patterns, as the 
latter time series had a similar time span, sampling and num¬ 
ber of data points. In the case of the correlation method, they 
are in a fairly good agreement (5.5 MEauh for the present time 
series versus 6.8 MEanh for the time series from observed pat¬ 
terns), whereas there is a certain discrepancy for the LPA method 
(9.9 MEarth and 15.7 MEanh for the simulated and observed RV 
time series, respectively). 


http://exoplanets.eu 

^ see Table 2 in Paper IV, case with no correction, total RV time se¬ 
ries, Set 1. 


18 












S.Borgniet et al.: Using the Sun to estimate Earth-like planets detection capabilities. 









Fig. 17. Detection limits versus stellar inclination for the correlation-based (four top panels) and the LPA (four bottom panels) methods. The detec¬ 
tion limits are displayed for different samplings and for different noise levels (black: no added noise, red: 0.01 m s“* noise, green: 0.05 m s“* noise 
and blue: 0.1 ms“* noise). Solid line: detection limits computed in case 2. Dashed line: detection limits computed in case 1. 
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Comparison between the two methods - The coiTelation- 
based method gives lower detection limits than the LPA method. 
This is also in agreement with our result s from Paper IV. Not e 
that we also compared the two methods in iMeunier et al.l (12012h . 
this time on actual stellar data. Out of a 10 target sample, we 
found in most cases that the LPA method gave lower detection 
limits. However in these cases, the targets were massive A-F 
Main-Sequence stars with medium to high v sin i (from 7 to > 
175 km s '), some of them being young stars. They would then 
have a very different activity pattern than the solar one. 

Impact of stellar inclination - For the LPA method, the 
decrease of the detection limit with the stellar inclination i in 
case 1 (“best-possible” detection limits) is consistent with the 
decrease of the activity-induced RV rms when going toward a 
smaller i. In case 2, we observe that the detection limits remain 
around 10 MEaith for i > 40° before increasing toward higher 
masses for smaller inclinations. This should mean that for 
i > 40°, the decrease of the planetary-induced RV amplitude for 
a decreasing inclination of the system is counterbalanced by the 
decrease in the same time of the activity-induced RV jitter. For 
smaller inclinations, the decrease of the activity RV jitter is less 
pronounced and hence it becomes more difficult to detect the 
planet. This is consistent with the fact that the behavior of the 
activity-induced RV jitter with respect to i is not sinusoidal. 

For the correlation method, the behavior of the detection lim¬ 
its is rather different. In case 1, the detection limits decrease with 
i until they reach a plateau at about 1.5 MEarth for i ~ 40°. When 
looking at case 2, surprisingly the detection limits begin to de¬ 
crease with a decreasing i, to reach a minimum and best value of 
2-2.5 MEarth for / = 40 - 50°. For smaller inclinations, the detec¬ 
tion limits increase with a decreasing i. Thus, the optimal config¬ 
uration for the detection of the planet is an orbital plane inclined 
of 45° with respect to the line of sight (in the most probable case 
of spin-orbit alignment). 

4.4.3. Impact of the parameters 

In this section, we study in more details the impact of the main 
parameters for the computation of the detection limits (i.e., the 
added RV noise level and the temporal sampling. To do so, we 
first study their influence on the detection limits computed above 
(Fig. [13. Then, we extend the study to a wider range of noise 
level or temporal sampling, but this time concentrating on one 
case (/ = 50°, case 2). This roughly corresponds to the inclination 
for which we get the better detection limits (in case 2). 

Temporal window - Our detection limit computation is based 
on the periodograms of the RV time series taken eight months 
a year, so as to better mimic actual observations. The peri¬ 
odograms will thus be modified compared to the periodograms 
of the total RV time series taken over the full simulation time 
span, due to aliasing. We compare the periodograms of the total 
RV time series (in the i - 90° case) before and after removing 4 
months a year in Fig. [18] Removing four months a year indeed 
injects a large amount of power in the RV periodogram around 
one year (~ 360-400 days) as well as at lower harmonics (180 
and 90 days). It will then impact our detection limit computation 
and deteriorate our detection limits, but it has to be taken into 
account to reproduce better the observations. 



1 10 100 1000 
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Fig. 18. Tot. RV periodograms for i = 90°, (1; 1) sampling. Black: 
data taken eight months a year. Red: full data. The 1% FAP are 
displayed (dashed lines), as well as the solar rotational period 
and its two first harmonics (blue), and the detection limit period 
(green). 

Temporal sampling - We found in Paper IV that when taking 
into account all RV components, the detection limits got worse 
for the (1;20) sampling (i.e., for a largely degraded sampling), 
while they did not vary significantly when going from a (1:1) to 
(1:4) or (1:8) samplings. Fig. [m shows that the LPA detection 
limits do not vary significantly when going from the (1:1) to the 
(1:20) sampling. As for the correlation detection limits, they do 
not vary significantly when going from the (1:1) to the (1:8) 
sampling, but get worse for the (1:20) case. 

We now focus on the i = 50° case to study the sampling im¬ 
pact in more details. We compute the detection limits for an ex¬ 
tended sampling range (from (1:1) to (1:70)) and display the re¬ 
sults in Fig. [T3 We find that both the correlation and the LPA 
detection limits are nearly independent from the temporal sam¬ 
pling up to the (1:10) sampling, and get worse for more degraded 
samplings (even if they show a large dispersion in the LPA case). 
The large dispersion at large temporal samplings probably re¬ 
flects the larger uncertainty on the detection limits for a smaller 
number of points. 

A(d(de(d noise - The detection limits computed with both the 
correlation and LPA methods are nearly independent from the 
added RV noise for a noise level up to 10 cm s ' (Fig.llTli. i.e. for 
the best RV accuracy expected on future spectrographs. This is 
in agreement with the results of Paper IV. We consider that the 
small variations seen between the detection limits computed for 
the different noise levels in the case of the correlation method 
are not significant. 

We study then the impact of the added noise for an extended 
range of noise levels, up to 1 ms' (i.e., the current RV accuracy 
reached on the better spectrographs such as HARPS), in the i = 
50° case (Fig.l20ll. The LPA detection limits are constant up to 1 
ms ', while the correlation detection limits remain constant up 
to 0.7 m s“' with a slightly larger value at 1 m s“'. 

5. Conclusion and perspectives 

We built a fully parametrized model of the activity pattern of a 
solar-like star, including the dark spots and the bright features 
(large faculae and network). The model includes about 30 pa- 


20 









































S.Borgniet et al.: Using the Sun to estimate Earth-like planets detection capabilities. 



Fig. 19. Detection limits vs. temporal sampling for i - 50°, no 
noise. Squares: LPA method; Stars: correlation method. 


8 


E 

c 

u 

v 

q5 

Q 


q\ -^^-.- 1 -^^^^^^^^-,- 1 -^^^^- 1 

0.0 0.2 0.4 0.6 0.8 1.0 

RV noise level (m/s) 

Fig. 20. Detection limits vs. RV added noise for i - 50°, (1:1) 
sampling. Squares: LPA method; Stars: correlation method. 


rameters that account for the different activity scales and the ac¬ 
tive structure behavior (most of them being well constrained by 
the litterature), and has a daily timescale over a complete activ¬ 
ity cycle. Using the same approach as in the previous papers, 
we deduced the corresponding RV and photometric time series, 
taking the inhibition of the convective blueshift into account in 
the case of the RV. The simulated activity pattern, as well as the 
time series, were compared to the work done in Papers 1,11 and 
III with data from solar observations. We found our model to be 
in remarkably good agreement with the previous data, thus as¬ 
sessing its validity. We then study the case of a solar-like star 
seen under different configurations so as to estimate the impact 
of the stellar inclination on the activity-induced jitter. Our results 
are the following: 

- For the stellar irradiance, the stellar inclination i has almost 
no effect on the amplitude of the long-term irradiance varia¬ 
tions. The decrease of the active structure filling factor with 
a decreasing i mostly counterbalances the enhanced contrast 
of the structures, in agreement with som e previous studies 
dKnaack et al.L 1200 ll I Shapiro et al.L l2014t) . On the contrary, 
the short-term jitter is reduced by a factor 6 when going from 
an edge-on to a nearly pole-on configuration. 


- For the RV variations, the peak-to-peak amplitude as well as 
the rms are reduced by a factor 8 to 10 on both the long-term 
and the short-term. For a nearly pole-on configuration, the 
remaining RV jitter is about 0.6 ms ' when considering the 
whole cycle, and is about 0.2 m s“* on shorter timescales. 

- When computed over the whole activity cycle, the activity- 
induced periodograms show power mainly in two period 
ranges: first at the stellar rotation period and its two or three 
first harmonics (induced by the crossing of active structures 
on the visible stellar hemisphere), and then at much longer 
periods, induced by more complex cycle effects. When going 
toward smaller inclinations, the power at the rotation period 
is reduced and nearly disappears for a pole-on configuration. 
The long-term cycle effects become preponderant. 

- In the case of a solar-like star, the convective component 
widely dominates the RV variations. The “photometric” con¬ 
tribution of spots and bright features account for about 35% 
of the RV jitter and decreases to about 15% for a nearly pole- 
on configuration. 

- For a solar-like star seen in a pole-on configuration, the 
photometric and RV (as well as presumably the chromo¬ 
spheric) activity-induced variations are most probably not 
distinguishable from those of a less active or almost inactive 
star. 

- Finally, when considering a realistic orbital configuration 
(f.e.,spin-orbit alignment), the optimal configuration for 
planet detection is a system inclined by about 45°. In this 
case, the lowest detection limits reach planetary masses of 
about 2MEarth at 480 days, without applying a correction to 
the RV signal. 

According to previous studies, solar- like stars show a grea t 
diversity of activity levels and properties (ISchroder et al.L 1^131) . 
The Sun itself is considered to be an average star, not partic¬ 
ularly active but not quiet either. Being fully parametrized and 
validated for the solar case, our activity model now allows us 
to explore a wide range of activity parameters and stellar prop¬ 
erties. In this paper, we have focused on solar-like activity and 
convection level stars; we expect stars with a lower activity level 
and/or stars with a lower convection level (such as K-type stars) 
to be much less affected and then to exhibit lower detection lim¬ 
its in the corresponding HZ. We also emphasize that activity 
correction methods (such as the Calcium index or the Ha line) 
are commonly applied to RV time series and already allow to 
reach detection limits at the MEauh level in the HZ (Paper IV). 
Optimized observation and re duction strategies are also promis¬ 
ing (such as averaging, see e.g. lDumusaue et ^1201 lah . Despite 
not being the focus of our paper, improving such methods and 
strategies, as well as looking for new ones, is an extremely im¬ 
portant question. Future instruments with very high RV accuracy 
may therefore be critical in implementing very efficient correc¬ 
tion tools to extract low mass planet signals. Our goal is to be 
able to test and even reproduce the photometric and RV vari¬ 
ations corresponding to each activity and stellar configuration. 
This should allow a better understanding of the magnetic stel¬ 
lar activity and open the way toward a better correction of the 
activity-induced stellar jitter. 

In Paper III, we derived the astrometric time series cor¬ 
responding to the observed solar activity pattern. Our activity 
model also allows us to produce the astrometric time series cor¬ 
responding to the simulated activity pattern (apart of the photo¬ 
metric and RV ones). We will present our results in astrometry 
in a separate paper. This is quite justified as upcoming instru¬ 
ments such as the Nearby Earth Astrometric Telescope (NEAT, 
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see e.g. iMalbet et al.L l2Q12h should allow for the first time to 
search for low-mass planets around nearhy FGK stars with as¬ 
trometry. 
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very useful comments on the manuscript. 


References 

Albrecht, S., Winn, J. N., Johnson, J. A., et al. 2012, ApJ, 757, 18 
Baliunas, S. L., Hartmann, L., Noyes, R. W., et al. 1983, ApJ, 275, 752 
Barnes, J. R., Jeffers, S. V., & Jones, H. R. A. 2011, MNRAS, 412, 1599 
Baumann, I. & Solanki, S. K. 2005, A&A, 443, 1061 
Berdyugina, S. V. & Usoskin, I. G. 2003, A&A, 405, 1121 
Cadavid, A. C., Lawrence, J. K., & Ruzmaikin, A. A. 1999, ApJ, 521, 844 
Cegla, H. M., Shelyag, S., Watson, C. A., & Mathioudakis, M. 2013, ApJ, 763, 
95 

Cegla, H. M., Watson, C. A., Shelyag, S., & Mathioudakis, M. 2014, ArXiv e- 
prints 

Chapman, G. A., Cookson, A. M., & Dobias, J. J. 1997, ApJ, 482, 541 
Chapman, G. A., Cookson, A. M., Dobias, J. J., & Walton, S. R. 2001, ApJ, 555, 
462 

Chapman, G. A., Dobias, J. J., & Arias, T. 2011, ApJ, 728, 150 
Claret, A. & Hauschildt, P. H. 2003, A&A, 412, 241 

Crouch, A. D., Charbonneau, R, Beaubien, G., & Paquin-Ricard, D. 2008, ApJ, 
677, 723 

Desort, M., Lagrange, A.-M., Galland, F., Udry, S., & Mayor, M. 2007, A&A, 
473, 983 

Dumusque, X. 2014, ApJ, 796, 133 

Dumusque, X.,. Boisse, I., & Santos, N. C. 2014, ArXiv e-prints 
Dumusque, X., Santos, N. C., Udry, S., Lovis, C., & Bontils, X. 2011a, A&A, 
527, A82 

Dumusque, X., Udry, S., Lovis, C., Santos, N. C., & Monteiro, M. J. P. F. G. 
2011b, A&A, 525, AMO 

Frohlich, C. & Lean, J. 1998, Geophys. Res. Lett., 25, 4377 
Fzresz, G., Epps, H., Barnes, S., et al. 2014, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 9147, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 9 
Galland, R, Lagrange, A.-M., Udry, S., et al. 2005, A&A, 443, 337 
Hall, J. C. & Lockwood, G. W. 2004, ApJ, 614, 942 
Harvey, K. L. & Zwaan, C. 1993, Sol. Phys., 148, 85 
Hatzes, A. P. 2002, Astronomische Nachrichten, 323, 392 
Hillenbrand, L., Isaacson, H., Marcy, G., et al. 2014, ArXiv e-prints 
Howard, R. E. 1991, Sol. Phys., 131, 239 
Howard, R. E. 1992, Sol. Phys., 137, 51 
Isaacson, H. & Eischer, D. 2010, ApJ, 725, 875 
Ivanov, E. V. 2007, Advances in Space Research, 40, 959 
Jeffers, S. V., Barnes, J. R., Jones, H. R. A., et al. 2014, MNRAS, 438, 2717 
Knaack, R., Fligge, M., Solanki, S. K., & Unruh, Y. C. 2001, A&A, 376, 1080 
Komm, R. W., Howard, R. F, & Harvey, J. W. 1993, Sol. Phys., 147, 207 
Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s 
grid. Kurucz CD-ROM No. 13. Cambridge, Mass.: Smithsonian 
Astrophysical Observatory, 1993., 13 
Lagrange, A.-M., Desort, M., & Meunier, N. 2010, A&A, 512, A38 
Lagrange, A.-M., Meunier, N., Desort, M., & Malbet, F. 2011, A&A, 528, L9 
Lanza, A. R, Bonomo, A. S., Moutou, C., et al. 2010, A&A, 520, A53 
Lanza, A. R, Pagano, L, Leto, G., et al. 2009, A&A, 493, 193 
Malbet, F., Goullioud, R., Lagage, P.-O., et al. 2012, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 8442, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 
Martinez Pillet, V., Moreno-Insertis, R, & Vazquez, M. 1993, A&A, 274, 521 
Megevand, D., Herreros, J.-M., Zerbi, R, et al. 2010, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 7735, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 
Meunier, N. 2005, A&A, 442, 693 

Meunier, N., Desort, M., & Lagrange, A.-M. 2010, A&A, 512, A39 
Meunier, N. & Lagrange, A.-M. 2013, A&A, 551, AlOl 
Meunier, N., Lagrange, A.-M., & De Bondt, K. 2012, A&A, 545, A87 
Noyes, R. W., Weiss, N. O., & Vaughan, A. H. 1984, ApJ, 287, 769 
Pepe, R, Mayor, M., Rupprecht, G., et al. 2002, The Messenger, 110, 9 
Rice, J. B. 2002, Asteonomische Nachrichten, 323, 220 

Saar, S. H. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 
294, Scientific Frontiers in Research on Extrasolar Planets, ed. D. Deming & 
S. Seager, 65-70 

Saar, S. H. 2009, in American Institute of Physics Conference Series, Vol. 1094, 
15th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, ed. 
E. Stempels, 152-161 


Saar, S. H. & Donahue, R. A. 1997, ApJ, 485, 319 

Santos, N. C., Bouchy, R, Mayor, M., et al. 2004, A&A, 426, L19 

Santos, N. C., Gomes da Silva, J., Lovis, C., & Melo, C. 2010, A&A, 511, A54 

Schrijver, C. J. 2001, ApJ, 547, 475 

Schroder, K.-P, Mittag, M., Hempelmann, A., Gonzalez-Perez, J. N., & Schmitt, 
J. H. M. M. 2013, A&A, 554, A50 

Shapiro, A. L, Solanki, S. K., Krivova, N. A., et al. 2014, ArXiv e-prints 
Triaud, A. H. M. J., Cameron, A. C., Queloz, D., et al. 2014, in lAU Symposium, 
Vol. 299, lAU Symposium, ed. M. Booth, B. C. Matthews, & J. R. Graham, 
399M00 

Ward, R 1966, ApJ, 145, 416 
Wright, J. T. 2005, PASP, 117, 657 


22 






